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An experiment has been conducted to acquire data for the validationof computational fluid dynamics codes used 
in the design of supersonic combustors. The flow in a supersonic combustor, consisting of a diverging duct with a 
single downstream-angled wall injector, is studied. Combustor entrance Mach number is 2 and enthalpy nominally 
corresponds to Mach 7 flight. The primary measurement technique is coherent anti-Stokes Raman spectroscopy, 
but surface pressures and temperatures have also been acquired. Modern design of experiment techniques have 
been used to maximize the quality of the data set (for the given level of effort) and to minimize systematic errors. 
Temperature maps are obtained at several p lanes in the flow for a case in which the combustor is piloted by injecting 
fuel upstream of the main injector and one case in w hich it is not piloted. Boundary conditions and uncertainties 
are characterized. 
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Nomenclature 

thermal conductivity of wall, W/mK 

number of samples 

number of parameters 

heat flux, W/m 2 

temperature, K 

time, s 

time of heater start, s 

time of fuel injection start, s 

position coordinates in right-handed system, m 

thermal diffusivity of wall, m 2 /s 

standard deviation 


Introduction 

C OMPUTATIONAL fluid dynamics (CFD) codes are exten- 
sively employed in the design of high-speed airbreathing en- 
gines. CFD analysis based on theReynolds-averagedNavier-Stokes 
equations uses models for the turbulent momentum, energy, and 
mass fluxes that employ many ad hoc assumptions and empirical 
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coefficients. Typically, these models cannot be applied with confi- 
dence to a class of flow for which they have not been validated. This 
experiment is one of a pair of experimentsconducted at NASA Lan- 
gley Research Center that have been adopted by NATO Research 
and Technology Organization working group 10 (scramjets) as a test 
case for their CFD developmentand validation activity. Much of the 
material reported herein, and associated CFD work, is reported in 
the literature, 1-4 as well as details of the other experiment, a study 
of a supersonic coaxial jet. 5 6 

The experiment was designed to provide a relatively simple case 
for CFD codes involving supersonic fuel injection, mixing, and com- 
bustion in a duct. The model geometry is simple, and large regions 
of subsonic/fecirculating flow are avoided, although the model con- 
tains a small step ahead of injection to assist with flameholding. 
The flow at the exit of the facility nozzle that provides the test gas 
(the supersonic heated air “simulant”) to the inlet of the model has 
been studied, both previously 7 and herein, including CFD calcu- 
lations, pitot surveys, one-dimensional analysis and coherent anti- 
Stokes Raman spectroscopy (CARS) measurements of temperature 
(see thefollowingdiscussion). Model wall temperature and pressure 
measurements were acquired. The model geometry was verified us- 
ing quality assurance techniques. The enthalpy of the test gas is 
nominally equivalent to Mach 7 flight. It was believed, on the ba- 
sis of calculations performed, 3 that this was hot enough to produce 
mixing-limited flow, that is, one for which chemical reaction of the 
injectedfuel proceedsat a much greaterto equilibriumrate than mix- 
ing. The advantage of a mixing-limited flow is that it is insensitive 
to chemical kinetics and, thus, provides a more unambiguous test of 
CFD turbulence models. However, it was later found that significant 
combustion of the fuel did not occur until a point well downstreamof 
injection, at which point a rapid rise in pressure and temperature was 
observed. The implication is that, between injection and this pointof 
“flameholding” chemical reaction greatly lagged mixing, whereas 
downstreamof this point the flow was indeed mixing limited. 

CARS is the primary experimental technique employed in this 
investigation. An introduction to CARS is given by Eckbreth, 8 and 
an application of CARS to supersonic combustors is given by Smith 
etal. 9 In the present application, molecular nitrogen is probed to ob- 
tain temperature. CARS, like many other optical techniques, does 
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not significantly alter the flow being studied. Intrusive probes, such 
as pitot, total temperature, hot-wire, and so forth, are not used due 
to access difficulty and high heat flux in the combustor, and because 
they may alter the flow. CARS has several ad vantages over other op- 
tical methods.lt is relatively mature and well understood. The signal 
is in the form of a coherent (laser) beam and is collected through rel- 
atively small windows. Incoherent (non-CARS ) interference is then 
rejected by spatial filtering. Consequently, signal-to-noiseratios are 
typically high when compared to techniques such as Rayleigh and 
Raman scattering. Furthermore, thermometry based on Rayleigh 
and Raman techniques often assumes the knowledge of pressure, 
which is not known in supersonic flow, 10 whereas CARS is rela- 
tively insensitive to modest pressure variations. Imaging techniques 
such as planar laser-induced fluorescence 1 1,12 or filtered Rayleigh 
scattering 13 offer instantaneous,quantitativemaps of temperaturein 
high-speed flows. However, a desirable aspect of CARS is that it is 
more suitable for multiparameter measurements suitable for com- 
piling statistics about the flow. These statistical correlations can be 
used to develop turbulencemodels for CFD. For example, the CARS 
system described in this paper is currently being modified to allow 
simultaneousmeasurementof temperatureand species mole fraction 
and, hence, to allow quantification of temperature- species correla- 
tions. CARS has been used to instantaneously measure temperature 
along a line. 14 However, due to the high-temperature, low-density 
conditions in the present experiment, we restricted measurements 
to a single point to maximize the signal intensity. 

Application of a complicated technique like CARS in high-speed 
engine environments is not routine. Because it is a pointwise (rather 
than planar) technique, building a “picture” of the internal tempera- 
ture field of the combustor requires hundreds of facility runs, which 
is expensive. These problems of experimental complexity and large 
data sets were mitigated by use of modern design of experiments 
(MDOE) techniques. Specifically, due to the complexity of the ex- 
periment there are likely to be significant uncontrolled variables, 
includingthose associated with facility operation, the model, instru- 
mentation, and soon. 15-17 Randomization techniques 18 are used to 
minimize systematic errors associated with these variables. Surface 
response methods 19-22 are used to represent the spatial dependence 
of temperature in each of a number of planes by a mathematical 
function of the spatial coordinates. These methods optimize the 
quantity of data that needs to be acquired to meet specific goals for 
the random error of temperature surfaces and provide a compact, 
functional form for these surfaces. Additional details of the appli- 
cation of MDOE to this experiment are described by Danehy et al. 2 

Flow Facility 

The experiment was conducted in NASA Langley’s direct- 
connect supersonic combustion test facility (DCSCTF; infor- 
mation available online at http://wte.larc.nasa.govyfacilities.cfm? 
id=2 [cited 29 September 2003]). Hot air simulant, known as vi- 
tiated air, is produced in the “heater” shown in Fig. 1. Molecular 
oxygen and air are premixed and H 2 is burned in this 02 -enriched 
air. Flow rates are selected so that the mass fraction of 0 2 is the 
same as that of standard air. The high-pressure vitiated air is accel- 
erated through a water-cooled convergent-divergent nozzle before 
entering the test model. The facility is designed to test the combus- 
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tor of a supersonic combustion ramjet (scramjet) engine by directly 
connecting the facility nozzle exit to the entrance of the combustor. 

The test conditions are nominally representative of Mach 7 
flight. Gas flow rates to the heater are 0.9153=0.008 kg/s for air, 
0.02843=0.0006 kg/s for H 2 , and 0.3003=0.005 kg/s for O 2 . The 
heater stagnation pressure is 0.765 3= 0.008 MPa. All uncertainties 
presentedin this paper are based on the 95% probability limits (1 .96 
times the standard deviation). The aforementioneduncertaintiesare 
due to the random run-to-run variations and do not include 3=3% 
uncertainty in the mass flow rate measurements. 

Heater and nozzle exit conditions are estimated from the flow 
rates, heater pressure, and nozzle minimum and exit areas using 
an inviscid one-dimensional analysis. The flow exiting the heater 
into the nozzle is assumed to be in thermodynamic equilibrium but 
has unknown enthalpy due to heat lost to the structure and cooling 
water. The enthalpy is guessed and the area of the sonic throat is 
computed by one-dimensionalanalysis, assuming isentropic flow in 
the nozzle. The enthalpy is then iterated until the computed area at 
the sonic throat equals the geometrical minimum area of the nozzle. 
Nozzle exit conditions are computed from the geometrical exit area. 
The composition at the nozzle throat and exit could be evaluated at 
frozen (at heater values) or equilibrium conditions. All significant 
minor species are included. Calculations assuming equilibrium and 
frozen composition differ in minor species concentration but not 
significantly in major species, temperature, or pressure. The nom- 
inal calculated conditions, and uncertainties due to mass flow rate 
measurement error and run-to-run variations in heater conditions, 
are heater stagnation temperature, 1 827 3= 75 K; exit temperature, 
1 187 3=60 K; exit pressure, 1003= 1.5 kPa; and exit Mach number, 
1 .989 3= 0.005 . Errors arising in the calculation due to the assump- 
tion of one-dimensional flow (the effects of nonuniform velocity, 
composition, temperature, etc.) are believed small. 

A study of the flow quality at the exit of the facility nozzle was 
previously conducted. 7 A pitot probe rake was employed to map the 
pitot pressure and the flowfield was visualized with the facility run- 
ning as a freejet (the model removed). Silane (SiH 4 ) was added to the 
heater H 2 and burned to form silica particles in the heater. The par- 
ticles were illuminated by a pulsed laser-sheet placed at the nozzle 
exit and imaged with a charge-coupleddevice camera. Results were 
compared to CFD calculations of the nozzle flow. Whereas the pitot 
pressure varied at the nozzle exit plane, the computed pitot pres- 
sure distribution agreed well with measurement. The flow appeared 
qualitatively more uniform in composition following a modification 
to the configuration of the facility’s ignitor. 

The test model is shown in Fig. 2; flow direction is from left 
to right. The model consists of two main sections of duct: a cop- 
per upstream section and a carbon- steel downstream section. Stain- 
less steel flanges and carbon gaskets separate the sections from 
each other and the nozzle. Proceeding from left to right, there is 
a constant-area segment, a small outward step at the top wall, and 
a second short constant-area segment followed by a segment with 
constant 3-deg divergence of the top wall. The span is constant at 
87.88 mm. Five small pilot fuel injector holes are located ahead 
of the step, and the main fuel injector is located just downstream 
of the start of the 3-deg divergence. The injection angle is 30 deg to 
the opposite wall. The injector nozzle is designed by the method of 
characteristicsto produceMach 2.5, one-dimensional flow at the in- 
jector exit. Molecular hydrogen injection is provided at a stagnation 
pressureof2. 12 3=0.07 MPa, a stagnation temperature of 302 3=4 K, 
and an equivalence ratio of 0.99 =t 0.04. On some runs, additional 
H 2 injection is provided by the pilot injectors at the same stagna- 
tion temperatureand a total equivalence ratio of 0. 148 3= 0.008. The 
pilots are turned on and off at the same time as the main fuel injector. 

The duct is uncooled; however, the wall thickness of the copper 
duct is greater than 32 mm and the carbon steel duct is 19 mm thick. 
Thus, given the good thermal conductivity of these materials, it is 
possible to operate the facility with the model fueled for run times 
in excess of 20 s (and unfueled for much greater times) without 
reaching excessive temperatures. With atmospheric-temperatuie air 
flowing in the model between runs, runs could be repeated every 
10-15 min. 
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Fig. 2 Test model: a) nozzle, copper and steel duct sections; b) detail 
in vicinity of fuel injector and pilots. Dimensions are in millimeters. 



B = Beam expander 
C = Combining mirror 
SL = Spherical lens 
S = Signal separator 
F = Filter 
P = Polarizer 
CL = Cylindrical lens 
E = Entrance aperture 
X = Splitter 


Periscope 


Fig. 3 CARS system: beam-combining optics; intersection and signal 
detection optics. 


The model is equipped with seven slots to allow the CARS beams 
to penetrate the duct, of which slots 1, 3, 5, 6, and 7, depicted in 
Fig. 2a, are used in this study. The slots are paired, one on each side 
of the duct. They are 4.8 mm wide and extend the full height of 
the duct. When not in use, the slots are plugged flush to the wall. 
Windows covering the slots are mounted at the end of short rectan- 
gular tubes at the Brewster angle to minimize reflections (Fig. 3). 
The window tubes are ventilated with a constant flow of electrically 
heated (~400 K) dry air to prevent condensation of water from the 
heater on the windows. The CARS interaction region can thus be 
translated the full span and height of the duct without damaging the 
windows. 

The model is instrumented with both pressure taps and wall tem- 
perature probes. Thirty-five static pressure taps are located in the 
copper duct, consisting of 0.80-mm-diam square-edgedholes. Taps 
are located on the bottom wall, at the centerline, and on the top 
wall at z = — 36.3 mm (z is measured from the horizontal center- 
line). Forty-nine static pressure taps are located in the steel duct, 
consisting of 1 .6-mm-diam square-edgedholes. Taps are located at 
top and bottom wall centerlines and sidewall midpoints. Pressures 
are measured with an accuracy of ±0.6 kPa by a Pressure Systems, 
Inc., electronically scanned pressure-measuring system. 


Six spring-loaded, bayonet, ribbon, K-type thermocouples from 
Nanmac Corp. are located in 6.35-mm-diam blind holes (from the 
outside) in the copper block. The bottom of each of these holes is 
square and the wall thickness (between hole base and duct interior 
surface) is 2.8 mm. A thermocouple junction is located at the bot- 
tom of each hole at the center. This choice of diameter and wall 
thickness is based on a multidimensional heat transfer analysis so 
that the thermocouple measures the unperturbed duct surface tem- 
perature. (While the temperature at the bottom of the hole is less 
than local duct surface temperature, the effect of the hole, which 
alters the conduction path, is to raise the local duct surface temper- 
ature. These effects offset.) Six “eroding” K-type thermocouples 
from Nanmac Corp., 7.9-mm diam and incorporatinga carbon-steel 
sheath, are located in the carbon-steelduct. The thermocouplejunc- 
tion for these probes is located flush with the duct flow surface. Over 
time, the junction was removed by the hot flow but was periodically 
regenerated by rubbing with sandpaper. 

CARS Technique 

Optical System 

The CARS system uses a Spectra-Physics DCR-4 unseeded, 
pulsed Nd:YAG laser, frequency doubled to 532 nm. The nomi- 
nal power is 550 mJ per pulse in the green, repetition rate is 1 0 Hz, 
and nominal linewidth is ~ 1 cm - 1 . A broadband dye laser was also 
used, consisting of an oscillator cavity formed between a total re- 
flecting mirror and a 30% reflector, and a single amplifier stage. 
Two identical Brewster- angle dye cells are used, one in the oscilla- 
tor and one in the amplifier, through which flows rhodamine 640 in 
methanol. The wavelength is centered between 605 and 606 nm, by 
adjusting the dye concentration, to match the Raman shift of nitro- 
gen. Figure 3 shows a schematic of the remaining components of 
the CARS system. The dye beam is expanded in a beam expander. 
Dye and 532-nm beams are combined with a dichroic mirror and 
are relayed via a periscope to a spherical (focusing) lens of focal 
length 0.4 1 m. All three beams cross at their focal points: the focal 
point of the dye beam is made coincident with the focal point of the 
532-nm beams by adjusting the beam expander. Beams are phase 
matched in a vertical planar BOXCARS configuration. 8 

At the lens, the full (to M0% of peak) diameter of the 532-nm 
beams is ~8.5 mm, and the beams are 1 8 mm apart. The full diam- 
eter of the dye beam is ~1 1 mm. At the focus, the diameters are, 
respectively, ~0.12 and ~0.15 mm, yielding a geometrical beam 
intersection length of 5.4 mm. The length of the measurement vol- 
ume was measured by translating the CARS measurement volume 
through a thin planar jet of nitrogen, surrounded by a coflowing 
jet of helium. The length over which there was any CARS signal 
is ~4.5 mm and the full- width half maximum is ~2.25 mm. The 
beam energy level per pulse at the focusing lens is ~85 mJ for each 
of the green beams and between 12 and 24 mJ for the dye. 

The beams (including the CARS signal beam) are relayed via a 
second spherical (collimating) lens and a second periscope back to 
the optical bench. The overlapping CARS and 532-nm beams are 
separated in a splitter. The splitter consists of two 100-mm-long 
dichroic mirrors that reflect ~99.5% of the incident CARS signal, 
while reflecting only ^20% of the 532-nm beams and transmitting 
the rest. The CARS and 532-nm beams enter the splitter at 45 deg to 
the mirrors and undergo six back-and-forthreflections. The beam is 
then directed through additional filters as needed to reject residual 
532-nm beams or for CARS signal attenuation. It enters a polarizer 
that allows only p-polarized light to pass (not a critical component 
because it is found that the signal is already well polarized). It 
is then focused by a pair of cylindrical lenses and enters a 1-m 
monochrometer with 1 200 groove/mm grating via an aperture. The 
cylindrical lenses are adjusted to produce a horizontal focus at the 
detector (maximizing image sharpness) and a loose vertical focus 
(to ensure the detector provides a linear response to signal 23 ). An 
optical splitter is located in front of the detector. 24 The splitter creates 
a secondary signal on the detector identical to the primary but offset 
and of 6% the intensity. When the intensity of the primary signal 
exceeds the dynamic range of the detector, the secondary signal is 
used for analysis. An EG&G PAR model 1420 intensified, linear, 
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self-scanned silicon photodiode array detector is mounted at the exit 
plane of the monochrometer. It consists of 1024 elements, 2.5 mm 
high by 25 /xm wide, and utilizes a fast (^millisecond decay rate) 
P-46 phosphor. The detector was scanned at 30 Hz, three times for 
each laser firing. The first scan was immediately after the laser pulse 
and acquired CARS signal, the second was used to clear the detector 
of any residual signal, and the third was used to determine the signal 
baseline. 

The two top prisms of the periscope are mounted on stepping- 
motor-driven vertical translation stages. The two bottom prisms and 
the vertical translation stages are mounted on similar horizontal 
stages. By translatingthe vertical and/or horizontal stages in tandem 
(maintaining alignment) the measurement volume is moved in the 
y and/or z direction. It is important that the beams are parallel to the 
direction of motion of the translation stage. During alignments, a 
reference, low-power helium- neon laser (He-Ne) beam is directed 
through the periscopes parallel to the direction of motion of both sets 
of stages and through the centerof the sphericallenses. The Nd: YAG 
and dye laser beams are periodically adjusted to ensure that they 
intersect the He-Ne beam at their focus. CARS data acquisition is 
under the control of a personal computer. Two types of acquisition 
are employed. In the first, data are acquired at a single point in 
space. In the second, data are acquired while either the vertical or 
the horizontal stages are in constant- velocity motion. 

CARS data are acquired in the supersonic combustor during mul- 
tiple sets of test runs. During a set of runs (which might last as long 
as 5 h), access to the model and optical system is prohibited. Test 
runs consisted of approximately 5 s during which the heater is op- 
erating but no fuel is injected in the model, followed by 1 1-20 s of 
fuel injection. CARS data are acquired over a period 2 s shorter than 
the period of fuel injection. Immediately after a run, 10 s of data is 
acquired with the dye laser beam blocked by a remotely operated 
flag. These “background” scans measure non-CARS interferences 
such as scattered laser light. 

Immediately before and after a set of test runs, “reference” CARS 
spectra are acquired in a low-speed jet of “PC-Duster” refrigerant 
gas. Because this jet contains no nitrogen, and is of constant CARS 
susceptibility, the spectra reflect the spectral variation in dye laser 
power. Additionally, spectra are acquired with room-temperatureair 
flowing in the duct. 

Large temperature variations occur in the test cell because it is 
continuously ventilated with external atmospheric air. Significant 
heat is also radiated from the model. Thermal expansion effects on 
the mechanical components of the optical system are not sufficient 
to require remote adjustments. (Adequate signal levels are typically 
maintained for several hours.) However, thermal expansions signif- 
icantly change dye laser power and center wavelength. Because it is 
possible for effects such as these to generate time- varying bias er- 
rors, MDOE quality assurancetactics were used. These tactics were 
designed to defend against such errors, to the extent allowed by con- 
straints in the operation of the apparatus, which will be discussed 
further. 

Data Reduction 

CARS data are analyzed on a workstation, separate from the data 
acquisition computer, after the data are acquired. Background scans 
are averaged and subtracted from data scans. Both primary and sec- 
ondary (producedby the splitter) CARS signal s are contained within 
the data scan. If the primary signal counts exceed the linear regime 
of the detector, the secondary is selected for analysis. Data scans 
are divided by the reference spectrum to remove the effect of the 
dye laser spectral power distribution and normalized to unit area 
(primary or secondary). Data are compared to a library of similarly 
normalized theoretical spectra to determine the temperature and ni- 
trogen concentration. The pixel location of the start of the theoretical 
spectra is allowed to float ± 1 0 pixels from a nominal position in 0.3- 
pixel increments. The combination of temperature, concentration, 
and pixel location that produces the least mean square deviation 
between theory and data is selected. 

The absolute concentration of nitrogen can be determined from 
spectral fitting if the optical system has been calibrated and the laser 


intensities are measured on a shot-by-shot basis. 25 However, for the 
sake of simplicity, this was not done in the current experiment. In- 
stead, the relative concentrationof nitrogen can be determined from 
spectral fitting if the CARS nonresonant susceptibility is known or 
can be reliably estimated. However, spatial variations in the compo- 
sition of the test gas were not quantified in the current experiment, 
so the local nonresonant susceptibility could not be computed, nor 
could a spatially averaged value of the nonresonant susceptibility 
be estimated that would be within a factor of 2 of the local nonreso- 
nant susceptibility. Therefore, neither absolute nor relative nitrogen 
concentration measurements could be determined in the present ex- 
periment. Nitrogen concentration was used as a fitting parameter 
in the data analysis simply to allow the spectral code to improve 
the quality of the fit, allowing the temperature to be determined 
accurately; the resulting values of nitrogen concentration were not 
quantitative and are not presented in this paper. 

Theoretical CARS spectra are generated using the program 
CARSFT. 26 The combustion gases are assumed to be a mixture 
of nitrogen and nonresonant buffer gas, both having nonresonant 
susceptibility of 8.5 x 1 0“ 17 m 3 /J. The nominal static pressure is as- 
sumed to be 1 atm (101,325 Pa). In fact, as will be seen later, the 
pressure in the duct at the CARS measurement planes is as much 
as 20 kPa above or below this value. The error introduced by this 
assumption is estimated by calculating theoretical spectra over a 
range of pressures and temperatures and fitting to them assuming 
1 atm pressure. The difference between fit and true temperatures (fit 
minus true) is ~25 K per 10 kPa difference from 1 atm for T s near 
2000 K. The exponential gap model for colli sional narrowing of the 
Raman line shape is used: a Voigt model was also tried and found 
to make little difference for the test conditions. A 532-nm laser 
linewidthof 1 cm -1 is assumed, though computed temperatures are 
not sensitive to the 532-nm laser linewidth for the test conditions. 
An experimentally determined instrument probe function is used. 
The width of this probe function was varied for best fit to a room- 
temperature CARS spectrum acquired at the beginning of each day 
of runs. Temperatures from 1 50 to 3000 K, in increments of 25 K, 
and 29 nonuniformly spaced concentrationsof nitrogen are used in 
the library. 

The peak of the reference CARS spectrum shifts significantly 
during a set of runs. This shift is largely due to a shift in the dye 
laser wavelength. Thermal effects on mechanical components of the 
optical system, beam steering by the hot gases in the duct, and even 
translation of the periscope stages all can change the position of 
the CARS beam focus at the entrance to the monochrometer and, 
hence, the position of the spectra on the detector. Two techniques 
are employed to correct for these problems. The first is to allow the 
pixel location of the start of the theoretical CARS spectra to float 
in the fitting process, as described. The second is to estimate the 
reference spectrum from the run data themselves. The technique 
is inspired by the observation that the directly measured dye laser 
spectra could be accurately fit to a Gaussian function. 

The procedure is illustrated in Fig. 4. The first step is to analyze 
the data from a given run using the reference spectrum obtained with 
the refrigerant gas at the end of a set of runs. For a particular laser 
pulse, spectrum A is the best- fit theoretical CARS spectrum, B is the 
measured CARS spectrum (prior to normalization by the reference 
spectrum), and C, the ratio of B to A, is an estimate of the refer- 
ence spectrum for this particular shot. All such estimated reference 
spectra for a given run are averaged and fit to a Gaussian function 
for amplitude, center, and width. The region of fit is indicated in the 
figure and excludes a region around the bandheadof nitrogen (which 
tends to be noisy) and the region to the left of pixel 275 (which is 
affected by masking from the optical splitter). As a practical matter, 
it is necessary to exclude from the average individual spectra for 
which good fits to theory could not be obtained. Acceptable aver- 
ages for fitting could be found for almost all runs. The Gaussian 
fit D is then used as the reference spectrum for normalization of 
all data for the given run. The figure shows data from a laser pulse 
after normalization, and the best-fit theoretical spectrum. This pro- 
cedure is typically iterated twice, with the latest best-fit theoretical 
spectrum used to determine an improved reference spectrum. 
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A = CARS theory 



Fig. 4 CARS data analysis. 

The techniques used for acquisition and analysis of CARS data 
in the supersonic combustor were tested in a “Hencken” adiabatic, 
flat-flame burner burning H 2 in air. Equivalence ratio was varied 
and 100 laser shot data runs were performed. The mean measured 
temperatures were found to be within ±50 K of temperatures cal- 
culated from measured flow rates, assuming equilibrium chemistry 
(including minor species). Data were also acquired in which 1) the 
total laser power was varied from 200 to 550 mJ and 2), through 
the use of different neutral density filters, the signal in the primary 
was saturated, forcing use of the secondary in analysis. No trends 
are found with either of these variables, indicating that the nitrogen 
spectrum is not saturated by high laser powers and that the splitter 
device works well. Standard deviations of temperature for a run are 
between 60 and 100 K. The trend with increasing laser power is 
to lower standard deviations owing to increased signal-to-noise ra- 
tios. Thus, the uncertainty (95% probability) of a single-shot CARS 
temperature measurement is ~ ±200 K, whereas uncertainty due 
to errors that are fixed over a run is ~ ±50 K. A previously noted 
discrepancy 1 in which the mean was found to be M50 K high at 
high equivalence ratio was subsequently shown to be an error due 
to a nonuniform flame induced by operation of the burner at 30 deg 
to the vertical (because of space constraints), and not a CARS error. 

Design of Experiment 

It is inevitable in a complex experiment such as this that there are 
sources of systematic variationthat are not controlled. These sources 
might be unidentified, or it may not be possible to correct for them. 
One example of a systematic variation that cannot be controlled 
in this experiment is the steady increase in wall temperature of 
the model that occurs during each run. For example, if this were 
to result in a steady increase in temperature throughout the flow, 
and if the data were all acquired in constant- velocity scans taken 
from left to right, then temperature maps would indicate the flow 
temperature in error (compared to the mean) high on the right and 
low on the left. Other examples include instrumentation drift, and 
instrument and tunnel operator fatigue and learning effects. It is 


possible to defend against these types of error by acquiring data 
in random order of spatial coordinate. This randomization has the 
effect of converting systematic errors into random errors. Unbiased 
estimates of mean temperature may then be obtained simply by 
averaging measurements at a point taken at different times. 

Practical constraints prevent the spatial order of data points from 
being fully randomized. Limited stepper motor speed causes motion 
from one point to the next to take significant time: if it is required 
to move between randomly selected points during the course of a 
limited-time facility run, the number of data points acquired would 
be considerably reduced. Changing from one plane to the next re- 
quires shutdown of the facility and hardware changes. Thus, all data 
at a given plane are acquired in a consecutive series of runs. During 
a given facility run, data are acquired either at fixed points or dur- 
ing continuous constant-velocitymotion of the CARS measurement 
volume, as described in the following. However, the order of these 
runs and directions of the individual motions are randomized. 

Data were acquiredover 201 facility runs over 10 testdays. Except 
for one day, when laser beams clipped the edge of the duct window 
slots due to thermal expansion of the duct during runs, the vast 
majority of the data were found acceptable and were analyzed. At 
each plane, data were acquired at six or seven fixed points near the 
horizontal centerline, 6 s with only the heater operating followed 
by 18 s with both heater and fuel injected. Data were also acquired 
during 1 6 s of horizontal motion of the translation stages at 5 mm/s, 
or from 9 s to 18 s of 5 or 6 mm/s vertical motion, during which 
time fuel was continuously injected. Two cases were considered. 
In the first “unpiloted” case the pilot injectors were not operated 
and data were acquired at planes 1, 3, 5, 6, and 7. In the second 
“piloted” case the pilot injectors were operated and data were only 
acquired at planes 1, 3, and 5 (due to facility availability constraints). 
Figure 5 shows a three-dimensional cutaway view of the duct with 
the locations of the data points overlaid. Flow direction is from top 
left to bottom right. 

Surface response techniques were employed to estimate mean 
temperature surfaces from the CARS temperature data. The data 
at each measurement plane are fitted to a cosine series bivariate 
function of order 5 (plane 1 only) or order 6, with the number 
of fit parameters respectively 21 and 28. Fits are performed us- 
ing commercial software (TableCURVE 3D, data available online 
at http://www.systat.com [cited 29 September 2003]). The number 



Fig. 5 Cutaway view of data acquisition points. 
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of data points acquired per plane, n, is from 2000 to 4000, and the 
standard deviation of the data from the fit at the various planes, a , is 
from 1% to 304 K. The uncertainty of the fitted functions, given by 
1 .96cr y/(p/n) (Ref. 19), is between 36 and 59 K depending on the 
plane. Here, p is the number of parametersof thefit. This uncertainty 
in the fits is actually an approximation representing the mean over 
all points to which the response model is fit. The true uncertainty 
depends on the specific location in a plane and is somewhat lower 
near the center of the plane and somewhat higher near the edge. 
Lower fit uncertainties could be achieved by increasingthe number 
of data points acquired. However, this is not warranted given the 
previously described uncertainty in the CARS mean temperatures 
of ±50 K due to fixed errors. 

Note that if a more conventional approach had been adopted in 
which data were obtained at a few discrete spatial points, the num- 
ber of samples required at each point to obtain an uncertainty of 
±50 K would have been approximately (1.96 x 304K/50K) 2 = 142 
samples/point. This number would allow only about (4000 samples/ 
plane)/(142 samples/point) =28 measurement points per plane. 
A functional representation of the surface temperature distribu- 
tion over the measurement plane, with given uncertainty, pro- 
vides better spatial coverage and a more compact representation 
than the 28 discrete points obtained with the same single-point 
uncertainty. 

In addition to the estimates of mean temperature surfaces just de- 
scribed, surface estimates of the standard deviation of temperature, 
or root mean square temperature deviation, were obtained. The dif- 
ference between the temperature and the mean obtained from the 
surface fit is obtained at each of the measurementpoints. These tem- 
perature deviations are squared and the resulting data are fit at each 
plane to a cosine series bivariate function of order 5. Standard devi- 
ation of these fits is between 78,000 and 1 90,000 K 2 and the fit un- 
certainty is between 15,000 and 32,000 K 2 . Typical values of uncer- 
tainty expressed as a percentageare ±25, 000/300 2 x 100= ±28%. 
It should be emphasized that this uncertainty is caused by insuffi- 
cient data points: as n increases, the uncertainty becomes small. 

Random deviations of temperature from the mean were traced to 
several sources. The deviation that we are interested in is inherent in 
the flow and is due to the fluid dynamical turbulent fluctuations. An- 
other source is random error inherent in the CARS instrument. The 
standard deviation of this error was measured to be ~ 100 K. The un- 
controlled systematic temperature variations, converted to random 
deviations by randomizing the order in which data were acquired at 
the various spatial locations, also contribute to the temperature fluc- 
tuations. The magnitude of these systematic variations is unknown 
but assumed to be small. Since these various sources of tempera- 
ture deviation are statistically independent, the mean square devia- 
tions are additive. Thus, in addition to the large random error, the 
experimental fitted mean squared deviation surface is everywhere 
~(100 K) 2 high. Given these high levels of uncertainty and error, 
the mean squared deviation surface fits are only qualitativelyuseful. 

Results 

Surface Pressure and Temperature 

Surface pressure and temperature data are presented for two typ- 
ical runs, one unpiloted and one piloted. These runs are chosen 
because the heater pressure and the gas flow rates to the heater, in- 
jector, and pilot are all very close to the respective averages over the 
total set of runs. 

Surface pressure distributions for the pressure taps on the top and 
bottom walls on the centerline, averaged over a 1-s interval at 16 s 
into the run, are shown in Fig. 6. Uncertainty in these measurements 
varies widely, depending on location in the duct, but are typically 
~ ±3 kPa due primarily to unsteadiness in the flow. Pressures vary 
widely in the upstream region due to the shock wave system created 
by the injectors and step. In the unpiloted case pressure generally 
falls, moving downstream due to divergenceof the duct, until 0.4 m, 
where it rises rapidly, peaking at ~0.7 m. Pressures for the piloted 
case are higher between the pilot (a = 0.074 m) and a* = 0.7 m due 
to combustion of gas from the pilot and main fuel injector in this 
upstream region. Downstream of 0.7 m the pressure drops smoothly 



Fig. 6 Surface-pressure distributions along centerline of top and bot- 
tom walls. 



Fig. 7 Surface-pressure time histories at selected taps on top wall for 
the unpiloted case. 


in both cases but is higher in the unpiloted case despite the lower 
total injected fuel rate. Differences between the two cases suggest 
that significant combustion of the fuel does not take place upstream 
of 0.4 m in the unpiloted case, but that the combustion of the fuel and 
air that are mixed then proceeds to completion by 0.7 m. It believed 
that fuel and air are not fully mixed at 0.7 m, so further mixing and 
combustion are expected to occur downstream. 

Surface pressure time histories for the unpiloted case, on the top 
wall, for several pressure taps between 0.4 and 0.7 m (in the re- 
gion of rapid pressure rise) are shown in Fig. 7. The heater flow 
is initiated at time t 0 = Is and fuel injection commences at time 
t j = 6.4 s. It is seen that, after the initiation of heater flow, there are 
large fluctuations in the pressure but no significant trends. Thus, the 
location at which combustion takes placeis not significantly affected 
by the increase in surface temperature that occurs during the course 
of a run. Much higher fluctuations in pressure are observed at taps 
downstream of the x = 0.464 m tap, suggesting that the combus- 
tion process is associated with a large increase in the level of local 
unsteadiness in the flow. The measured pressure is filtered by the 
response of the tubes between the surface and the transducer, and 
the response function of the transducer. Because the system time 
response is much longer than characteristic turbulence time scales 
of the flow, these data underestimate the unsteadiness. 

Because the duct is uncooled, surface temperatures vary greatly 
during the course of the run, as seen in Fig. 8 for selected points in 
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Fig. 8 Wall temperature history, unpiloted, top wall: data and fits. 



x, m 


In the piloted case, q \ is large on the top wall, downstream of the 
main injector (jc = 0.166 m), indicating the presence of combustion 
products near the wall . In the unpilotedcase, q x is significant (though 
not so large as the piloted case) ahead of 0.7 m (where the pressure 
rise begins), consistent with incomplete combustion near the wall. 

CARS Temperatures 

Figures 10 and 1 1 are three-dimensional cutaway views of the 
duct showing contour plots of the fitted temperature functions and 



400 x 

300 


Fig. 10 Cutaw ay view of duct showing contour plots of mean temper- 
ature for unpiloted and piloted cases. 


Fig. 9 Surface heat flux distributions: top and bottom walls, 
centerline. 

the unpiloted case. In the copper section, temperature is typically 
~320 K at the start of fuel injection but rises to as high as ~500 K 
in the unpiloted case and ~590 K in the piloted case. In the carbon- 
steel section, temperature is typically ~320 K at the start and as 
high as ~890 K at the end. These variations are least-squares fit 
to the solution for wall temperature T w of a semi-infinite body at 
initially uniform temperature T, and subject to steps in surface heat 
flux from zero to q 0 at time t 0 and from q 0 to q 0 +q j at t\ (Ref. 27): 



The fit parameters are T ’• , q 0 , and q x and the fit is conducted out to 
/ = 1 1 s. As seen in Fig. 8, fits often diverge from the data beyond 
1 1 s, indicating that heat flux continues to change slowly during the 
run. The material property yj(x s / k s is taken to be 36.7 kW • K/m 2 s 1 ^ 2 
for the copper duct and 12 kW • K/nrs 1 / 2 for the steel duct. Heat 
fluxes are presented in Fig. 9, both top and bottom walls, on the 
centerline. Heat flux varies from point to point in the range 0.3 to 
1 .0 MW/m 2 before injection (q 0 ), and from 0.7 to 1 .8 MW/m 2 with 
injection ( q 0 + q [ ), where q x is the increase in heat flux. At a given 
point, q 0 varies 10 to 30% from run to run, probably reflecting 
variation in the initial temperature of the wall between runs. At 
downstream locations, is large due to injection and combustion. 



Piloted 


Unpiloted 


Fig. 1 1 Cutaway view of duct showing contour plots of rms tempera- 
ture deviation for unpiloted and piloted cases. 
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the square root of the fitted mean-square deviation functions (i.e., 
the rms deviation). Also shown are two of the four duct walls. Re- 
call that the main fuel injector is on the top wall between planes 
1 and 3, and that the pilot injectors are on the top wall upstream 
of plane 1. An idea of the scale may be found by referring back to 
the dimensioned drawings in Fig. 2. Please note that the tempera- 
ture measurements, although essentially instantaneous, are widely 
spaced in time relative to any characteristic time of the flow, and so 
consecutive measurements are uncorrelated. Thus, it was not possi- 
ble to extract from these measurements such quantitiesas turbulence 
time scales. 

Consider first the fitted mean temperature surfaces (Fig. 10). At 
plane 1, in the unpiloted case, the temperature is fairly uniform, 
between 1030 and 1250 K. The mean temperature of this plane is 
1 162 K. This mean compares favorably with the value computed as- 
sumingone-dimensionalflowfromtheheater,whichis 1 187±60K. 
In the piloted case, the fitted temperature is also reasonably uni- 
form, except that the temperature drops slightly close to the top 
wall, downstream of where the pilot fuel is injected. There is no 
indication of pilot fuel combustion in this plane. At plane 3, in the 
unpiloted case, there is a region with temperature as low as ~250 K 
at the center, which is the injected fuel plume. There is no evidence 
of combustion of the injected fuel. In the piloted case, there is a 
band of hot pilot fuel combustion products close to the top wall. 
The center of the main fuel plume may be seen as a cooler region 
with temperature as low as ~650 K, greater than in the piloted 
case, suggesting some combustion. Plane 5 is similar to plane 3. 
Temperatures have risen near the center of the fuel plume, to as 
low as ~550 K in the unpiloted case and ~1250 K in the piloted 
case. The height of the region of hot combustion products near the 
top wall in the piloted case has increased. At planes 6 and 7, in 
the unpiloted case, temperatures have risen abruptly compared to 
plane 5, suggesting nearly complete combustion (i.e., H 2 and0 2 do 
not coexist at any given point). In the regions in these planes with 
the highest temperatures (~2300 K), the fuel and air are probably 
mixed to nearly stoichiometric proportions. The cooler region near 
the center (M500 K) is fuel rich. No significant quantities of in- 
jected fuel have penetrated to the sidewalls (M 200 to 1 300 K). Data 
are not acquired at planes 6 and 7 in the piloted case due to lack 
of time. 

Figure 12 shows temperature measurements from a single ver- 
tical scan near the horizontal centerline for plane 5 in the un- 
piloted case. The temperature surface fit along this spatial line is 
also shown, with error bars, as are the mean temperatures from 
several fixed-point runs, also with error bars. The error bars of the 
surface fit are based on all the data in the plane and so include the 
effects of run-to-run variations. The error bars of the fixed points 
are based on the standard deviations of the given run, containing 
about 150 individual measurements. In neither case is uncertainty 
due to fixed errors included, which was given as ±50 K. It may be 
seen that the fixed-point mean temperatures agree with the surface 



Fig. 12 Temperature data at z- 3.5 mm, plane 5, unpiloted case. 


fit, providing confidence in the approach adopted and suggesting 
that run-to-run variation is not significant. This figure also illus- 
trates graphically the relative levels of uncertainty and the high lev- 
els of fluctuations in the data. Note that the fluctuations either side 
of the temperature surface fit indicate real temperature fluctuations 
in the flow, because fluctuation levels are larger than the instru- 
ment precision error of ~ 100 K (standard deviation). The presence 
of high temperature fluctuations should not be surprising since, at 
many points in the turbulent combusting flow, there will at cer- 
tain times be hot combustion products and at other times cooler air 
or fuel. 

Consider now the rms temperature deviation surfaces (Fig. 11). 
All the surfaces have irregular oscillations of amplitude MOO K 
due to random error. However, there are significant features present. 
At plane 1, for the unpiloted case, the rms deviation surface varies 
between ~ 100 and ~230K. The standard deviations obtained from 
measurements at fixed points varied from 140 to 220 K, significantly 
greater than that inherent in the instrument (MOO K). Thus, there 
are significant real temperature fluctuations in the flow exiting the 
heater, with standarddeviationM 00 K or greater. The rms deviation 
level in the piloted case is comparable, except that it rises close to 
the top surface, downstream of the pilots. 

Several regions of large rms temperature deviation levels may be 
observed. One region of high rms deviation level, up to M40 K, 
lies between the center of the fuel plume and the top wall in plane 5 
of the unpiloted case. A second region, with rms deviation level up 
to MOO K, lies near the fuel plume in plane 3 of the piloted case. 
In both cases, these regions appear to be ones where combustion 
is just beginning or is incomplete. This suggests that intermittency 
in degree of combustion is the main cause of the fluctuations; that 
is, in some lumps of fluid passing through the measurement plane 
combustion has taken place and in others it has not. The result for 
the unpiloted case is particularly interesting because there is little 
indication in the mean temperature plot of combustion having taken 
place, but solid evidence in the rms plot of intermittent combustion. 
One further region of large rms deviation level, with peak rms de- 
viation level of ~525 K, lies on the right side near the top wall, in 
plane 5, in the piloted case. The mean temperature is M200 K at 
this location, much higher than for the other locations discussed, 
so the rms deviation level represented as a fraction of the mean is 
much lower than in these locations. It is probable that the temper- 
ature surface fit has greater error at the edge of the data plane. The 
high standard deviation at this location may simply be explained by 
greater error in the mean temperature surface. 

The fitted mean temperature and rms temperature deviation sur- 
faces, and the surface-pressureand heat flux data, provide a consis- 
tent description of the flowfield. In the unpiloted case, only small 
and intermittent combustion of the injected fuel is observed ahead 
of x =0.4 m. Downstream of x =0.7 m, combustion appears nearly 
complete. In the piloted case, combustion of the pilot fuel appears to 
take place between x = 0. 122 (at the step) and 0.274 m. There also 
appears to be significant combustion of the main injected fuel by 
x = 0.274 m. Theseresults are not consistent with CFD calculations 
performed prior to the commencement of the experimental work, 3 
which predicted combustion in the vicinity of injection in both un- 
piloted and piloted cases. Consequently, this experiment provides a 
test case for CFD that is more challengingthan anticipated. Accurate 
calculation will require accurate modeling of both the chemical ki- 
netics and the turbulence- chemistry interactions,as well as accurate 
modeling of the turbulent mixing. 

Conclusions 

CARS thermometry and MDOE techniques have been success- 
fully applied to a simple supersonic combustor. Functions that pro- 
vide an estimate of the mean temperature and rms temperature de- 
viation surfaces at various planes in the flow have been obtained. 
Additional data include surface pressure and temperature distribu- 
tions. Boundary conditionsand uncertaintiesare well characterized, 
and the data are useful for the validation of CFD codes used in the 
design of supersonic combustors. Accurate calculation of the ag- 
gregate data set will require accurate modeling of both the chemical 
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kinetics and the turbulence-chemistry interactions, as well as accu- 
rate modeling of the turbulent mixing. 
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